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In order to clarify thermodynamic stability of Mg-based long-period stacking ordered (LPSO) 
structure, we systematically study energetic preference for alloys on multiple stacking with different 
composition for random mixing of constituent elements, Mg, Y, and Zn based on special quasirandom 
structure (SQS). Through calculation of formation free energy of SQS, Mg-Y-Zn alloy exhibits phase 
separation into Mg- and Y-Zn rich phase, which is consistent with previous theoretical studies. Bulk 
modulus of SQSs for multiple compositions, stacking sequences, and atomic configurations ranges 
around 35 GPa, i.e., they do not show significant dependence of Mg concentration, which therefore 
means that the effects of phonon do not play significant role on LPSO phase stability. Introducing 
stacking fault to hep stacking gains ’’negative” energy, which indicates profound relationship between 
introducing stacking faults and the formation of long-period stacking ordering. 


1. INTRODUCTION 

For its remarkable high tensile strength and ductil¬ 
ity a, Mg-based long-period stacking ordered (LPSO) 
structures are considered as light-weight structural al¬ 
loy for next generation. In order to clarify relationship 
between thier formation process and resultant proper¬ 
ties in terms of application for structural materials, con¬ 
siderable number of experimental as well as theoretical 
studies have been carried out. Previous theoretical stud¬ 
ies mainly address formation process and thermodynamic 
stability of Mg-based LPSO alloy, including (i) the ten¬ 
dency of phase separation confirmed by Cluster variation 
method [2|, (ii) in-plane ordering of clusters consisted 
of Y, Zn substitutional atoms [3|, and (iii) systematic 
understanding of energetic stability with respect to va¬ 
riety of substitutional atoms into Mg-based alloys Q. 
Although these previous theoretical works partly clarify 
themodynamic stability of LPSO phases, they did not 
sufficiently discuss about (i) relative stability in terms 
of disordered phases or (ii) effect of lattice vibration 
on stability of LPSO. The former one is considered es¬ 
sentially important, since well-established ordering en¬ 
ergy, determining the thermodynamic stability of ordered 
phase with respect to temperature, is typically reffered 
to the difference in mixing energy between ordered (here, 
LPSO) and disordered phases. For the latter one, vibra¬ 
tional effects, their significant role on phase stability has 
been amply demonstrated for several binary alloys, such 
as Ref. [9[. Fig. [1] is a schematic illustration of stack¬ 
ing sequence of 18R LPSO structure, which includes Mg 
hep stacking and Y, Zn substitutional atom concentrated 
phase including LI 2 cluster. To systematically evaluate 
thermodynamic stability of such structures, we need to 
consider energetics of multiple structures including dif¬ 
ferent stacking sequence. However, most calculations 
confine the structural model, so calculation that doesn’t 
confine experimentally reported structures is highly re¬ 
quired. Furthermore, it is fundamentally important to 
assess thermodynamic stability of LPSO phase by com¬ 



FIG. 1. Schematic illustration of stacking sequence of 18R 
LPSO structure. Stacking sequence is for hexagonal closed 
packed (hep), which are denoted by A,B and C. ”S.F.” means 
stacking fault. 


peting disordered phase. In this study, based on DFT 
calculation, we systematically study thermodynamic sta¬ 
bility of LPSO structure in terms of disordered phases. 


2. METHODOLOGY 

To comprehensively address energetics of multicompo¬ 
nent system with various composition, we need to cal¬ 
culate tremendous number of structures of the system 
and computational cost therefore become large. In this 
study, we employ special quasirandom structure (SQS) 
to assess thermodynamic stability of LPSO phase com¬ 
peting with disordered phase. SQS is a special micro- 
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scopic state whose multiple correlation functions are nu¬ 
merically identical to those in perfect random structure, 
which therefore provides physical properties for perfect 
random alloy [lOj. As following procedure, we calcu¬ 
lated correlation function of SQS of ternary system. 

Let us consider the system with N lattice points for 
number of components, R. Gi is a variable which specifies 
the occupation of lattice point i and a = ( 01 , 02 , • • •, gn) 
can specify any atomic arrangement and we can represent 
structures in following equation, which is called correla¬ 
tion function; 


= PdA°i)Pd 2 {vj) ■ ' ' Pd n (<J k ). (1) 

Here, p(cq) is complete ortho normal basis func¬ 
tion at lattice point i and it is obtained by applying 
Gram-schmidt technique to the linearly independent 
polynominal set ( 1 , cq, of, • • •, cr^ _1 ). a denotes a 
cluster included in the structure, d n is the index of basis 
function p(cq) and M is the set of index. In the case of 
ternary alloy system, the occupation of lattice point by 
each element, Mg, Y, and Zn, is defined bycr = +l,—1,0 
respectively, leading to the basis function of eq. (| 2 ]) & 


Po = 1 ,pl = y|cr, P 2 = -\/2(l - ^CT 2 ). (2) 

By averaging over equivalent clusters in lattice, 
correlation function of SQS in Mg-Y-Zn ternary system 
with compositions of each element, XMg, xy and xzn, 
, is represented by eq. (J3]) , 


= 2 2 • 3 ri 2 m ( XMg -x Y ) n m (l-3x Zn ) m , (3) 

where n and m denote dimension of clusters and the num¬ 
ber of p 2 respectively. The simulation was performed so 
that correlation functions come closer to ideal SQS cor¬ 
relation functions up to 4-th nearest neighbor pair (i.e. 
6 pair clusters on hep) finded by eq. (J3]) . We optimized 
correration functions for each clusters by performing nu¬ 
merical simulation 8 ] and constructed SQSs. To eval¬ 
uate the accuracy of obtained correlation functions of 
SQS, we compared the simulated values of correlation 
functions with standard deviation of them in configura¬ 
tional space. Then, most of the errors of correlation func¬ 
tion are small enough compared with standard deviation. 
Based on these obtained structures, we calculate forma- 

d 2 E 

tion free energy, Ff orm , and bulk modulus, B = —V——r 

ov z 

to evaluate thermodynamic stability of Mg-Y-Zn system 
that has multiple stacking sequence with various com¬ 
position for random mixing. We constructed structures 
with various composition on mutiple stacking sequence, 
whose composition is shown in Figure. [ 2 ] The detail of 


calculation condition of structure is shown in Appendix. 
Fform is denoted by eq. and we define U by eq. ©, 

Fform = U ~TS (4) 


U = E(Mg x Y y Zn z ) - xE{ Mg) - yE(Y) - zE{ Zn), (5) 

where E{^) is total energy of 7 , and we evaluate config¬ 
urational entropy, 5, based on Bragg-Williams approxi¬ 
mation. In order to evaluate bulk modulus, we calculate 
total energies of structures, whose volume is expanded at 
a rate of ±3%, ± 6 %, ±9%. By using B based on Debye- 
Griineisen approximation, Debye temperature, ©d, can 
be described as 0 


©d = 41.63 



M 


(6) 


where Vo is atomic volume and M denotes atomic mass. 
Based on empirical Debye model, F v ib can be estimated 
by eq. 0 ; 


F v ib = jk© D + kT{ 31n(l - exp(^) - D(^-)}, (7) 


where D(x) is Debye function. 

We estimated the effect of phonon on Mg-based alloy 
by calculating and the vibrational free energy of 

formation with respect to the pure constituents can be 
estimated by eq.(| 8 ]) [ 20 ]; 


ajttco = F vib (T) - ][>q ib (T), 


(8) 


where i denotes element and Xi is a composition of i. 

In order to evaluate the effect of stacking faults on sta¬ 
bility, we constructed a SQS which has stacking faults on 
hep (stacking sequence is ” ABABCABABCABABC”). 
Hereinafter, we call this structure ”mixed”. The ef¬ 
fect of stacking difference is quantified by interfacial 

-L'mixed Efcc,hcp 


energy, which is defined by Li.p. = 


A 


where F? m ixed is total energy of structure including stack¬ 
ing faults and F?f cc ,hcp is for structures on fee and hep 
stacking. A denotes area of interface, which is 1.4 xlO 2 
A 2 in this calculation. Additionally, to address stabil¬ 
ity of ordered phase competeing with disordered phase, 
we also calculated formation energy of structures in¬ 
cluding Y-Zn LI 2 cluster and estimated ordering tem¬ 
perature. We employ first-principles calculations using 
a DFT code, the Vienna Ab-initio Simulation Package 
(VASP) [ 12 ] [lH, to obtain the total energies for struc¬ 
tures of Mg-Y-Zn alloys. The calculation of total en¬ 
ergy is carried out for the structures in Table. IA11 All¬ 
electron Kohn-Sham equations are solved by employing 
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the projector augmented-wave (PAW) method [14] [15|- 
We select generalized-gradient approximation of Perdew- 
Burke-Ernzerhof (GGA-PBE) [16[ form to the exchange- 
correlation functional. The plane-wave cutoff energy is 
set at 350 eV throughout the present calculations. Bril- 
louin zone sampling is performed on the basis of the 
Monkhorst Pack scheme [lg. K-point mesh is set 4x4x1 
for structure, C3, and 4x4x4 for others and smearing pa¬ 
rameter is 0.15 eV Hil. 



FIG. 2. Composition of structural model. Further informa¬ 
tion are denoted in Appendix. 


3. RESULTS AND DISCUSSION 

First, we evaluated formation free energy of struc¬ 
tures with multiple compositions on hep and fee stack¬ 
ing, which is shown in Fig. [3l At T = 500K, Ff or m 
of structures on hep is negative at all through compo¬ 
sition and possesses two extreme values. Through calcu¬ 
lation of formation free energy of SQSs, Mg-Y-Zn alloy 
exhibits phase separation into Mg- and Y-Zn- rich phase 
and this result is consistent with previous research by 
Iikubo et al. Q, which suggests the validity of this sim¬ 
ulation based on SQSs. Using these optimized SQSs, we 
proceed our discussion to calculate bulk modulus, order¬ 
ing temperature, and the effect of stacking difference on 
phase stability. 

Next, in order to evaluate the effect of lattice vibration 
on phase stability of Mg-Y-Zn system, we calculated bulk 
modulus, which is shown in Figure. |U Dashed line repre¬ 
sents linear averaged bulk modulus, H ave . Bulk modulus 
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FIG. 3. (left) Formation free energy of SQSs with multiple 
compositions on hep and fee stacking with respect to energies 
of structures on hep. Curves are guide for the eyes, which 
shows tendency of phase separation, (right) Formation energy 
of structures including LI 2 cluster (F4, G2) and SQSs on hep 
(FI) and fee (F2, Gl) stacking. The ratio of concentration 
of substitutional atom is denoted in each figure. Structure 
indices are denoted in Table. lATl 


for Al, El HI, FI (on hep stacking), E2 (on fee stack¬ 
ing) and F4 (ordered structure) are about 35 GPa. The 
values of bulk modulus are smaller than F> ave through 
all compositions. Calculated AF^ m at T = 1000K at 
(xMg,^y,^Zn)=(87.5,6.25,6.25) was found to be negligi¬ 
bly small within the calculated accuracy of B. Compared 
with avove results of Ff orm , F v ib does not have signifi¬ 
cant influence on Ff orm of Mg-Y-Zn alloys with multiple 
stackings and compositions. Additionally, bulk modulus 
does not show significant dependence of Mg concentra¬ 
tion, which therefore means that the effects of phonon do 
not play essential role on LPSO phase stability. In this 
study, since effect of optical mode for multicomponent 
system is not considered, further study including opti¬ 
cal mode is needed to quantitively clarify the effect of 
phonon for stability of Mg-Y-Zn alloys. Hereinafter, in 
evaluateing effects of ordering and stacking sequence dif¬ 
ference on stability, we only consider the configurational 
effect. 

Then, we examine the effect of ordering on stability. 
We constructed structure including LI 2 cluster, F4 and 
G2, whose atomic arrangements are shown in Fig. [5] For¬ 
mation energies of F4 and G2 are lower than SQS on hep 
(FI) and fee (G1,F2) stacking (right side of Fig. [3]), which 
shows this system is stabilized by the effect of ordering. 
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xMg (at.%) 

FIG. 4. Bulk Modulus of SQSs and ordered structure. Dahsed 
line shows B ave = xBMg+yBY + zBzn , where Ba is bulk mod¬ 
ulus of A and x, y, z denotes composition of each element. 
Colors correspond to the ratio of concentration of substitu¬ 
tional atom. 
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FIG. 5. Schematic illustration of the atomic arrangements of 
structures including Y-Zn LI 2 clusters. Red and blue spheres 
represent Y and Zn atom respectively. F4 : Meta-stable. G2 
: Stable. 


As we estimated ordering temperature, T c , from F2 to 
F4 and from G1 to G2 was about 1050 K and 930 K 
respectively. This results shows that transition tempera¬ 
ture depends on the arrangement of LI 2 cluster and T c of 
G2 is close to experimental results by Okuda et.al. [211 


taking that transition temperature can be overestimated 
based on Bragg-Williams approximation into considera¬ 
tion. Moreover, this result indicates that Mg-Y-Zn sys¬ 
tem form ordered phase up to melting point: Mg-Y-Zn 
alloy can be regarded as inter metallic compounds. 

Finaly, we address the effect of stacking sequence on 
stability using interfacial energy. Fig. [6]is interfacial ener- 
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FIG. 6. The effect of stacking difference on interfacial energy 
of Mg-Y-Zn alloy on various stacking sequence 


gies of structures on different stacking (FI, F2, F3) with 
respect to the energy of structure on fee stacking, F2. As 
shown in Fig. [6l interfacial energy of disordered phase 
that introduces stacking fault is lower than the energy 
for linear average of hep and fee stacking. This result in¬ 
dicate that when stacking fault is introduced to Mg-rich 
hep alloys to form ”Ll 2 -like” ordering, corresponding in¬ 
terface between original hep and formed fee region gains 
” negative” energy, which is contrary to the conventional 
tendency that interface gains positive energy. This spe¬ 
cific characteristics for interface energy in Mg-based al¬ 
ios would thus be one of the fundamental prerequisite to 
accelaration of forming LPSO phases, which should be 
further investigated in the future work. 


4. CONCLUSION 

In this study, we calculate preference of energetics of 
Mg-Y-Zn ternary alloy system in terms of disordered 
phase stability. Through calculation of formation free en¬ 
ergy of SQSs, Mg-Y-Zn alloy exhibits phase separation 
into Mg- and Y-Zn- rich phase. Bulk modulus for SQSs 
and ordered structure ranges about 35 GPa and they 
do not show significant dependence of Mg concentration, 
which therefore means that the effect of phonon does 
not play essential role on LPSO phase stability within 
acoustic mode. Order-disorder transition temperature is 
estimated about 930 K and this results is colse to the 
melting point that is experimentaly reported and this 
suggests that Mg-Y-Zn LPSO alloy can be regarded as 
inter metallic compounds. The effect of stacking faults 








5 


stabilized the LPSO phase and this results indicate that 
there remains profound relationship between introducing 
stacking faults and the formation of long-priod stacking 
ordering. 
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Appendix 


In this study, we calculated free energy and bulk mod¬ 
ulus of structures shown in following data. 


TABLE Al. Composition, number of atoms, and stacking sequence of calculated structures 


index Mg (at.%) Y:Zn number of atoms stacking sequence model calculated bulk modulus (GPa) 


Al 



96 

hep 

pure Mg 

35.2 

A2 



96 

fee 

pure Mg 


B1 



96 

hep 

pure Y 

50.4 

B2 



96 

fee 

pure Y 


Cl 



96 

hep 

pure Zn 

34.9 

C2 



96 

fee 

pure Zn 


D1 

93.75 

1:1 

96 

hep 

SQS 


D2 

93.75 

1:1 

96 

fee 

SQS 


El 

87.5 

1:1 

96 

hep 

SQS 

35.5 

E2 

87.5 

1:1 

96 

fee 

SQS 

35.4 

FI 

85.4 

4:3 

96 

hep 

SQS 

35.5 

F2 

85.4 

4:3 

96 

fee 

SQS 


F3 

85.4 

4:3 

240 

mixed 

SQS 


F4 

85.4 

4:3 

96 

fee 

Order (LI 2 ) 

35.9 

G1 

80.6 

4:3 

72 

fee 

SQS 


G2 

80.6 

4:3 

72 

fee 

Order (LI 2 ) 


HI 

75 

1:1 

96 

hep 

SQS 

35.6 

H2 

75 

1:1 

96 

fee 

SQS 


I 

50 

1:1 

96 

hep 

SQS 




